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Stephen  D.    Lowe'1' 

Lawrence  Radiation  Laboratory- 
University  of  California 
Berkeley,    California 

October  30,    1963 

ABSTRACT 

A  full  solution  for  the  fractional  release  of  a  tracer  element 
through  a  moving  boundary  has  been  obtained  for  the  cases  of  a  slab 
and  a  sphere.      In  both  cases,    the  initial  distribution  of  the  tracer  is 
assumed  to  be  uniform  throughout  the  body.      For  a  constant  rate  of 
boundary  motion,    the  full   solution  is  valid  for  all  times  over  the   range 
from  zero  to  complete  evaporation  of  the  body,    and  this   solution  may 
be  applied  to  all  finite   rates  of  boundary  motion. 

By  applying  certain  limitations  to  the  time  range  and  to  the  total 
amount  of  boundary  motion,  the  full  solutions  are  reduced  to  approxi- 
mate forms  that  are  presented  graphically  as  a  family  of  curves. 

For  a  more  limited  range  of  application,    i.    e.  ,    for  small  values 
of  the  quantity    —  ( — J   (— —J     '      (<  0.2),     where    b   is  the   rate  of  boundary 
motion,    t   is  the  time,     D   is  the  diffusion  constant,    and   a    is  a  charac- 
teristic dimension  of  the  body   (half -width  for  the   slab  and  radius  for 
the   sphere),    simplified  expressions  are  developed  for  the  fractional 
release  of  the  tracer,    f,    which  may  be  used  to  more  accurately  deter- 
mine the  diffusion  coefficient,    D. 

These   simplified  expressions  are 

f=   -4-   (^)l/2+   i(^)    for  the   slab, 


and 


1/2    V   2/ 

'  a 

6        /Dt\l/2        3   /bt 


,  6       /Dt\l/2        3  /bt\  ,        .,  , 

f  =   7—   {  — )    '       +   —  (- — )  for  the   sphere. 

1/2    \ZJ  2   \a  /  ^ 

it     '  xa 


TT    '  a 
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I.     INTRODUCTION 
A.     Basis  for  the  Problem 

During  his  investigation  of  the  diffusion  of  xenon  in  uranium 

monocarbide,    Shaked  noted  the  slope  of  his  plot  of  the  fractional  re- 

/      1/2 
lease,     f   vs    (t)    '         increased  noticeably  during  a  high-temperature 

1 

anneal.       He  postulated  that  evaporation  of  the  uranium  monocarbide 

was  responsible  for  the  increase. 

1/2 
A  similar  increase  in  the   slope  of  the     f   vs    (t)    '         plot  had  been 

reported  previously  for  experiments  in  the  temperature  range,     1600 
to  2200°C. 2 

With  the  increasing  interest  in  high-temperature  applications 
for  uranium  fuels,    in  thermionics,    for  example,    it  is  considered  of 
interest  to  investigate  the  release  of  a  fission  product  from  the  fuel 
element  due  to  the  additional  process  of  evaporation  occurring  simul- 
taneously with  the  process  of  diffusion. 

B.  Discussion  of  the  Problem 

Several  investigators  have  reported  diffusion  constants  for 

3 
xenon-133  and  krypton-85  in  uranium  dioxide.       A  lesser  number 

have   reported  results  for  the  diffusion  constant  of  xenon-133  in 

1 
uranium  monocarbide.        In  the    "conventional"    analysis,    one  assumes 

a  uniform  initial  concentration  of  the  tracer  element  in  the  geometri- 
cal body  whose  boundaries  are  fixed  with  time.      Evaporation  of  the 
body  is    discussed  in  a  general  manner,    but  the  effects  are  minimized 

4-11  1.    3 

expe  rimentally . 

C.  Statement  of  the  Problem 

This  report  is  an  investigation  of  the  effect  of  evaporation  upon 
the  fractional  release  of  a  radioactive  tracer  element  from  a  fuel 
body.      As   such,    the  decay  of  the  tracer  is  neglected  for  the  purposes 
of  this  report    (i.  e.  ,    X  =  0),     although  this  effect  can  and  should  be 
taken  into  account  for  experiments  involving  long  postirradiation 
anneal     time. 

An  expression  for  the  combined  fractional  release  due  to  both 
diffusion  and  evaporation   (i.  e.  ,    diffusion  through  a  moving  boundary) 


-2- 


is  obtained.      This  result  is  compared  to  the  fractional  release  derived 
for  the  case  of    "pure  diffusion"   (i.  e.  ,    diffusion  through  a  stationary- 
boundary). 

The  analysis  is  performed  for  the  case  of  the  finite   slab  and  the 
case  of  the  finite   sphere.      It  is  assumed  that  the  initial  concentration 
of  the  tracer  element  in  the  body  is  uniform  and  that  the  concentration 
of  the  tracer  element  at  the  boundary   (moving  or  stationary)    of  the 
body  vanishes. 


II.      THE  SLAB  PROBLEM 
A.      Concentration  of  the  Tracer 
One   seeks  a  solution  for  the  diffusion  equation 


acC^^m**)         _(a-^t)<x<a-^  (ii-i) 


Jt 


such  that  C(X,0)=C, 


(11-la) 


and  CL±(Q-bt)lt']   =  0  ~for     0<t<^/bJ 


(II -lb) 


where    bis  the  rate  of  movement  of  the  boundary  that  may  be  obtained 

from  kinetic  theory  if  one  assumes  evaporation  in  a  vacuum  of  a 

4 
material  of  known  vapor  pressure.        The  assumption  of  a  constant 

evaporation  rate  ignores  the  possible  effect  upon  vapor  pressure,    and 
hence  upon  the  evaporation  rate,  if  there  is  a  change  of  chemical  com- 
position of  the  solid  surface  accompanying  evaporation.     At  time 

t  =  a/b,    the   sample  will  have  completely  sublimed, 

5 
Chambre  has  shown     that  the  solution  to  Eq.  ( II - 1)  that  satisfies 

an  initial  even  function    f(£)     and  Eq.  (II- lb)  may  be  expressed  as 


CfctW 


2(77- Dt/a 


A 


-f- 


&iTjf(s)**p  [t&^SCxt^cOJj] 


where  for  the  case  under   study 

and    SU,t)5)=  expl-^f]^  exp[-C*zg 


(H-2) 


(II-2a) 
(H-3) 


-4 


Making  the   substitutions  into  Eq.  (II -2),    one  obtains  the   result- 
ing expression  for    C(x,  t) : 


cfr«-^feej*Jie*p[-#]4 


-a. 


+M**r  f-V+pp^W^ 


(ii -4) 


m-i      i 


a- 


HDt 


Performing  the  integration  over    £   (see  Appendix  D),    one  may  show 
Eq.  (II -4)    to  be 


Cm-^-Urfj^-erFfefr] 


+ 


c 


-|(-0-exp{-K¥-^)]]HP(^)[erft(t(;iat) 


_^vr  a(2Y)-l)-(x+2M)  1 
€n       2(Dt)fc  J 


(II -5) 


+  ^|(-/f  exp{-[^^-^)J[  U(-^ a(^/>fr-^^ 


2  Cot)' 


ert     2CDty^       Jr. 


This  may  be  rewritten  as 


+ 


c 


(II -6) 


2©)' 


Making  the  definitions 


£- 


-  5t 


(II-7) 


-5- 


and 


where  0  <  6  T  <  -1; 

and  substituting  into  Eq„  (II -6),    one  obtains 

-     Co 


(H-8) 


(II-9) 


C 


*  §-|M"exp{-ryi[».(^t)-^i}  (erf  W^'V 

7LV1*-  j 

By  using  Eq.  (C-4)  and  setting    B=  0   (i.  e.  ,    no  evaporation  or  boundary 
motion),    Eq.  (11-9}  becomes 


er 


zz''*<      erfc  zz''*- 


(11-10) 


zz7*-  zz'z     J 


This  may  be   rewritten,    after  combining  terms,    as 

CfV)  =  C-C„fC-^[^c%^i+eH£%^].      <n-n: 

Equation  (11-11)  may  be  recognized  as  one  equivalent  to  Eq„  (A-6), 
which  shows  that  the  expression  for  the  tracer  concentration  does   re- 
duce to  that  for  the   stationary  boundary  case  when  no  evaporation  is 
considered. 
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B.      Total  Fractional  Release 

The  total  release  of  the  tracer  element  due  to  the  combined  ef- 
fects of  diffusion  and  evaporation  may  be  obtained  by  integrating  the 
concentration  at  any  time    (t)    across  the  volume  of  the   slab,    dividing 
the  result  by  the  total  initial  amount  of  tracer  present,    and  subtract- 
ing this  quotient  from   1.      Thus, 

£=  a.-Qft)  =  i -as?),  (u-iz, 

+  Qo  '  (So    ' 


where 


6}0  =   Q^(2a-Co)  (II-12a) 


and 


Qfr)-  c(\z)j(%)  ^  zjcMdH) 

Therefore,  (i-az) 

_r=   /_  zj0  ccva,r)d(%) 


(II -12b) 


-  J-  J ^|f  )n€xp{-^  K'-^)+x4]]  U^^')H^pt) 


£TV* 


Performing  the  integration  and  combining,    one  has 

2T^       f  in-") 

For    r  =  0(t  =  0),     the  fractional  release,    given  by  Eq.  (II- 14),    is 
zero,    as  is  to  be  expected;  and  for    (3t  =  1   (or  complete  evaporation), 
the  expression  becomes 

■f  =  1  +  Z"z  (Wc  ^k ~  Cerfc  ^  ] 

•     ^  .      r         -,vr»-l 

4 


^^T-W"'{^H*-«-'"^^ 


/2. 


which  is  equal  to   1,    after  cancellation  of  equivalent  terms,    and  cor 
responds  to  the  total  release  of  the  tracer. 

Setting    (3  =  0   (no  evaporation),    Eq.  (11-14)    becomes 


(11-15) 


This  may  be   reduced  to 


which  may  be  recognized  as  equivalent  to  Eq.  (A-8),    the  fractional  re- 
lease under  conditions  of  a  stationary  boundary. 

C.      Approximate  Expression  for  the  Fractional  Release 

For  the  usual  values  of   t    encountered  in  fission  gas -release 
experiments    (r  <  0.01),    the  expression  obtained  for  the  fractional  re- 
lease  in  Eq.  (11-14)  may  be  immediately  reduced  to 

-f  =  fSZ  +  Z'^itsrf^^ 

"■"''  ZT''X  (11-17) 

If  the  further  limitation  is  made  tha,t    (3t<  0.5,    Eq.  (11-17)  may  be 
reduced  further  to 


■fs/ST+'cVrft^ 


-9- 
of  which  the    n  -  1    term  is  of  the  greatest  significance,    as  follows: 

By  using  the  relations   (C-13)    through   (C-18),    Eq.  (11-19)  be 


come  s 


^+  Z^[,erfcf  +g(,r(^r>er{c  ^ } 


(11-20) 


St1/2                  6t1/2 
For  large    ILL- (i.  e.  ,    E-L- >2.5),     Eq.  ( 11-20)  reduces  to    f  =  (3t, 

which  indicates  that  the  release  of  the  tracer  is  due   solely  to  evapora- 
tion. 

For    (3  =  0,    Eq.  (11-20)  becomes 

f_      2,      r'/a,  (II-21) 

the  approximate  expression  for  fractional  release  from  a  slab  with 
stationary  boundaries. 

Figures  II.  1  and  II. 2  show  the  fractional  release  from  a  slab   vs 

a  /  "y  A  /  "J 

[time]    '      (t         )    for  various  rates  of  evaporation  ((3).     As  a  reference, 
the  stationary -boundary  plot  is  given,    and  in  each  figure  it  is  the 
lower  line.      Figure  II.  1  is  for  the  longer  periods  of  time,    whereas 
Fig.  II. 2  represents  shorter  periods  of  time. 

Figures  II. 3  and  II. 4  demonstrate  the  contribution  of  the  individ- 
ual terms  that  make  up  the  expression    [Eq.  (11-20)]    for    f.      Figure 
II. 3  represents  an  evaporation  rate  of   (3  =  10    for  times  less  than 
t  =  0.01,    and  Fig.  II. 4  covers  the  time  spanless  than   t  =  0.0001    for  a 
higher  evaporation  rate  of   (3  =  50. 
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(1) 

p  =       0 

(2) 
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(3) 

p  =      10 

(4) 

p  =     20 

(5) 
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(7) 
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no  evaporation) 


1/2 
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Fig.      II.  1.        Fractional  release 
of  evaporation. 
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for  various  rates 
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Fig.     II.  2.        Fractional  release    vs 
of  evaporation. 


1/2 


for  various  rates 
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f=PT+Zr1/2ierfc£ll/2  +  T1/2     £    (4)m  (V/2Km+l  erfc  SI* 

2  m=l  \2  /  2 


/2 


—  Release — no  evaporation  (p  =  0) 

-  -  Combined  curve  for    (3=10 


(4)    Tl/2(4)2(|Ll/2Ji3erfc|Ll/2 
,,.       1/2,..3/bt1/2"\3- 


(1)  Pt 

(2)  2r1/2ierfc^l/2 


.4       r     Sri/2 
i    eric  *— - 


MU.)19!) 


Fig.      II.  3.        Components  of  fractional-release  curve, 


13- 


o    l/2  •       r     (3T1/2         1/2     f-    ,  ,.mf  StVA™    m+1        .     Qt1/^ 
=  pr +  2T1/^ierfc  ^j-         +t    '         2^    (4)      I    ,  x  erfc 

2  m=l  \Z  /  2 


Release — no  evaporation 
Combined  curve  for    3  =  50 


(DPT  (3)    r1/2(4)(|Il/2)i2erfcfl/2 

(2)    Zr^ierfclL172       (4)    x1/2(4)2('f  l/2)V erfc  f U* 


0.012 


0.006 


0.008 


1/2 


Fig.      II.  4.        Components  of  fractional-release  curve. 
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III.      THE  SPHERE  PROBLEM 
A.      Concentration  of  the  Tracer 
A  solution  for  the  diffusion  equation  is  sought.      Let 

such  that  C(r,  0)=    Q  (III -la) 

and  C(^-tt,t)=    0  "for  0<t<^,  (Ilt-lb) 

where    b,    again,    is  the   rate  of  boundary  movement.      At  time    t  =  a/b, 
the   sample  will  have  been  completely  sublimed. 

Equation  (III - 1 )  may  be  reduced  to  that  for  the   slab  problem  by 
making  the   substitution 

UL(r>t)  =  rC(rjt)j  (in-2) 

subject  to  the  conditions 

u.(r,  o)  =  r  C0  (in -2*) 

and  U.^-ttJi):B    0  for  0<t<C^"    •  (III -2b) 

Since    C(r,    t)    is  finite  at  interior  points,    it  is  necessary  to  define 

u.(  0,-t)  -   0  •  (III"2c) 

Chambre      has  shown  the   solution  of  Eq.  (III-l),    which  satisfies  an 
initial  even  function    f(£,)    and  Eq.  (Ill- lb),    to  be 


(III-3) 


n-i   -&, 


-15- 


where,    for  this  problem 


#3)  -  C 


(III-3a) 


and     SfrltJ5)-eV[-^]-€xP[-^] 


(III -4) 


By  making  the   substitutions  into  Eq.  (III-3),    the  resulting  expression 
for    C(r,   t)    is 


+ 


(HI-5) 


-etp(-^^t)W- 


By  performing  the  integration  over    |    (see  Appendix  E),    Eq.  (III-5) 
may  be   shown  to  be 

+  ertc       a(Dt)"*       J 
^    cerfc2CDt)7-J 


-^(#[^c 


'4  r  fl-r     _  : prr^  ( 


(III-6) 
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Using  the  definitions  (II  -7)  and  ( II -8),    then  substituting  into  (III-6), 

one  obtains  . 

_  Co      r„~£.  l-^/cL        —r.  l  +  r/a 


Wfi-c-ffal^zv* 


-  er-fc 


ZX 


M 


+  er-fc  ZrrCt^hhlliz 

zz'1^ 

t  ey^c  2n()-i3z)- 1+1-/^1 
ZZ'1 


+ 


(in-7) 


'/i* 


Setting    (3  =  0   (i.  e.  ,    no  evaporation  or  boundary  movement)    in 
Eq.  ( III  -7),    it  follows  that 


2r 


■o 

Zr 


>?=/    L  ZZ  ZZ  h- 

-  c-^zuz  [^"fe^y*  -  ^^  ^y*- 


(III-8) 


■17- 


The      ierfc     terms  reduce  to  zero,    and  after  some  recombination, 
Eq.  (III-8)  becomes 


C(M  =  CD-^2[erfc£^-erfc&^ 


(III -9) 


This  may  be   recognized  as  equivalent  to  Eq„  (B-6),    which  shows  the 
expression  for  the  tracer  concentration  to  reduce  to  that  for  the 
stationary -boundary  case  when  no  evaporation  is  considered. 

B.      Total  Fractional  Release 

The  total  release  of  the  tracer  element  due  to  the  combined  ef- 
fects of  diffusion  and  evaporation  may  be  obtained  by  integrating  the 
concentration  at  any  time    (t)    throughout  the  volume  of  the   sphere, 
dividing  this   result  by  the  total  initial  amount  of  tracer  present,    and 
subtracting  this  quotient  from  1.      Hence, 


where 


and 


Thus, 


C-  Qo-Q(^)  .   /_  act) 


(111-10) 


(III-lOa) 


.   Q(Z)  =  \U~&(rLfC(%Z)d?k).  (in-iob) 

C  =  I  _  ^l7rl?C(r/aiz)d(r/«) 

T  W^  Co 


+ 


18- 


+  erfc^('^l-l+r^ld(^) 


\  2.  C  /7- 

_  ierfc  2dt^zI±rk.\<jQrL) 

=  /  -  3j?3w)  4lcS[^  ^-^c  ^]^) 


flSg^-r^^H^'1 


+erfc^^^?^)    (m-11) 


+  e  ^"c  M^llllhi  }  4(»L) 


+■  3  r  ,4|c£)  [*Hi  4j$t  -  ^rfc  4ggfe-  ]  <tf  %) 


19- 


Integrating  and  combining  the  terms  in  (III- 11),    one  finds 

zz'         J 


YA--1 


+  L™erfc&jW-£z)-l 


(111-12) 


m, 


y/t-  Z> 


yn~-3 
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For    T  =  0(t  =  0),    the  fractional  release,    as  given  by  Eq.  (Ill- 12),    is 
zero,    as  expected;  and  for    (3r  =  1   (or  complete  evaporation),    the  ex- 
pression becomes 

{  =  3  -  3  +1  -  <*  ?  (£**&  ar*  ~  C*erk  «W 
-,2Z3/Z(  i3erfc  jfe*.  ~  i3erfc  ^) 


which  is  equal  to   1,    after  cancellation  of  equivalent  terms,    and  cor 
responds  to  total  release  of  the  tracer. 

Setting    (3  =  0   (no  evaporation),    Eq„  ( III  - 12)  becomes 


f^fllW^  +  t'er^ 


(111-13) 


OO 


This  may  be   reduced  to 

(III  -14) 


f  =  6  TV*  [»"*  *  2 i^6rfc  ^  ]  ~  3  T , 
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which  may  be  recognized  as  equivalent  to  Eq„  (B-8),    the  fractional  re 
lease  under  conditions  of  a  stationary  boundary. 

C.      Approximate  Expression  for  the  Fractional  Release 

For  the  usual  values  of   r   encountered  in  fission-gas -release 
experiments   (t  <  0.01),    the  expression  obtained  for  the  fractional  re- 
lease in  Eq.  (Ill -12)  may  be  immediately  reduced  to 

+  C  Z  0 -pi)  ft rfc  |^  -  1 2  Z %  i3e rfc gfa 


m'z  (■  __  /2T'/a  (111-15) 

Yri  =  Z  r  2  X1'2- 

+  exp  i  nftp+iX1?*)]  terfcB»±!X!^ikl  J 

rn-3  Zz<u 

By  making  the  further  limitation  that    (3t  <   0.5,    one  obtains 

-IZZ^t3erfc-t^T^ 
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(111-16) 
n-i  /2  L 


of  which  the    n=  1    term  is  of  the  greatest  significance,    as  follows 


-  C  7(l-(3z)Z(-irfaz''f-zC>»erfi  =££ 


(111-17) 


By  using  the  relations    (C-13)    through   (C-19),    Eq.  ( III  - 17)  becomes 

f  =  3pz(l-pc)+$zf-  3T  ~  3Z(l-(3z) 

+  &T/z0-pt)ier(c  gf\  IZT0-f3T)izer-fc  4^ 
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+•  3Zk  {l-pc)i  W(fT  r'e  rfc  ^ 


(111-18) 


'/z\m  arr'k 


+  C.ZZ,<"»iX4r(&T  "r%rtc-& 


The  last  three  summations  will  be  small  compared  to  the  first  terms; 
therefore,     f    will  be  approximated  by 


^=3^(l-^)^f-3Z-^0-/3l)  +IZZ0-/3Z)  terfcg 


+  3Z,k0-^)[2ierfcM,\i  (V^^V^  JT 


^c  nn;i       '    ^  £. 


(111-19) 


R    d/2                    r    i/Z 
For  large     BL- (i.  e.  ,    BL- >  2.5),    Eq.  (Ill  -19)  will   reduce  to 

f-  Sp^(l-f3z)-T^f~-^Z-5Z(l'^), 

which  indicates  that  the  release  of  the  tracer  is  due  primarily  to 
evaporation. 

For    R  =  0,    Eq.  ( III  -  19)  becomes 

-f=4LZ*-3Z<  (ni-20) 

the  approximate  expression  for  fractional  release  from  a  sphere  whose 

boundaries  are   stationary. 

1/2 
If  a  further  restriction  is  made,    such  that    r    '      <0.0l   (or 

t  <  0.0001),    then  those  terms  in  Eq.  ( III  - 19 )  having    t   as  a  coefficient 

(excepting    Rt    terms)    may  be  neglected,    and  one  obtains 
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-f  =   5pZ  (1-fSz)  +<jet)  (m.21) 

f  3 T  * (I-/3Z)  [2 i€r&  ^  V  t  (iTqtjMrfcqp 


St*/2 

For  large    *- ,     Eq.  (111-21)  reduces  to 

-f=3p<r(l-(3?)+(/3T)3>  I 

which  indicates  the   release  of  the  tracer  to  be  due  solely  to  evapora- 
tion. 

For    (3  =  0,    Eq.  ( III  -21)  becomes 

-f=  £k  Z'/Z  >  (ni-22, 

the  more  familiar  expression  for  fractional  release  from  a  sphere 
with  stationary  boundaries. 

Figures  III.l  and  III. 2  display  the  fractional  release  from  a 

l/2        l/2 
sphere    vs    [time]    '      (t         )    for  various  rates  of  evaporation   ((3).      As 

a  reference,    the   release  for  the  case  of  a  stationary  boundary   ((3  =  0) 

is  given  for  each  figure;  see  line    (1).      Figure  III.l   (note  the  dotted 

line)    in  addition,    shows  the  effect  of  neglecting  the    3t    term  in 

Eq.  ( III  -20)    for  longer  periods  of  time.      Figure  III. 2   represents  the 

release  of  the  tracer  element  for  shorter  time  periods. 
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Fig.      III.  1.       Fractional  release    vs    t 
of  evaporation. 
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for  various  rates 
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(  1)  (3   =        0  (no  evaporation) 

(2)  (3  =     20 

(3)  p  =   100 

(4)  (3  =  200 


(5)  p  =     500 

(6)  (3  =  1000 

(7)  p  =  2000 
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Fig.     III.  2.       Fractional  release    vs    t   '    '    for  various  rates 
of  evaporation. 
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IV.      DISCUSSION 

A.      General  Results 

For  all  values  of   (3>0,    the  fractional  release  with  time  is  in- 
creased over  the  corresponding  value  for    f    when    (3  =  0.      For  any  value 
of  time    (or    t),    f    will  increase,    for  both  the   slab  and  sphere  cases, 
as    3    increases.      With  smaller  values  of   6,    say    3<  Z0,    the  curve  of 

0.0 
1/2 


f   vs    t  very  closely  approximates  a  straight  line  for    t    '     <0.01, 


which  indicates  that    f    is  proportional    (or  very  nearly  so)    to    t 

l/2 
The   same   statement  may  be  made  for    t         <  0.1,    provided    6<  2. 

B.      Initial  Release 

This  theory  does  not  account  for  the  large  early  release   rate  in 

postanneal  experiments  noted  by  several  observers.    '         In  these 

1/2 
cases,    the   slope  of  the    f   vs    t  curve  decreased  with  increasing 

time  until  a  region  of  more  or  less  constant  slope  was  attained. 

C.     Spherical  Case --Stationary  Boundary 

Figure   III.  1  illustrates  the  effect  upon    f    of  neglecting  the   second 

1/2 
term  of  Eq.  (111-20).      For  values  of    t         <  0.01,    the  difference  in 

values  of    f     resulting  from  lack  of  consideration  of  the   second  term  of 

1/2 
Eq.  ( III  -20)  is  negligible.      For  larger  values  of   t         ,    the  fractional 

release  predicted  by  Eq.  ( III  -22)  becomes  increasingly  greater  than 

that  predicted  by  Eq.  (111-20),    and  differ  s  by  0. 03    (~    10%)    at    t1/^2  =  0.1. 

Since  both  Eq.  ( III  -2  0)  and  Eq.  ( III  -22 )  appear  to  be  plotted  as   straight 

lines,    an  error  of  17%  in  the  value  of    D    results  for  the  above  value  of 

t    from  neglecting  this  term.      It  is  noted  that  Eq.  ( III  -22 )  closely 

represents  the  release  curve   resulting  from  boundary  motion  for 

which    (3  =  2. 

D.      Spherical  Case --Moving  Boundary 

1/2 
For  values  of   t         <  0.01,    Eq.  ( III  -21)  may  be  used  in  preference 

to  Eq.  ( III  - 19)  for  all  values  of    (3t.      If    (3t   (remember    (3t    independent 

of    D    and  is  bounded  between    0    and    1)    is  approx<0.2,    then  the  term 

3 
((3t)      may  also  be  neglected. 
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1       l/z 

For  large  values  of    -=- (3t         ,    the  fractional  release  is  approxi- 
mated by     3(3t(  1  -  0t)  +  (Bt)3>     which  is  a  cubic  equation  in  time    (or    r). 
This  may  occur  at  short  times  for  very  rapid  boundary  motion   (0),    or 
for  very  long  periods  of  time    (t).      Neither  of  these  cases,    however, 

is  useful  in  making  a  determination  of   D. 

1        1/2 
For  small  values  of    —  0t         ,    pertaining  both  to  relatively  short 

periods  of  time  and  small  rates  of  boundary  motion,    the  fractional  re  - 

1/2 
lease  is  directly  proportional   (or  very  nearly  so)    to    t  and  there  - 

1/2 
fore  to    t         ,      It  should  be  possible,    therefore,    to  develop  a  simple  ex- 
pression that  may  be  used  to  determine    D   with  greater  accuracy  than 

gTl/2       /aTi/Z  \ 

would  the  use  of  Eq.  (111-22).      For  small    *- — ,    I <  0.2)  ;  hence, 

Eq.  (111-21)  may  be  approximated  by 

With  Eq.  (IV- 1),    one  can  predict  the   release  curve  to  be  para- 
bolic in  shape  when    f    is  plotted   vs    time  to  the  one -half  power;    the 
curvature  depends  upon   b.       The  curve  consists  of  a  straight -line  com- 
ponent and  a  parabolic  component,  which  is  the  evaporation  correction. 

Figure  IV.  1  is  the  graphical  presentation  of  Eq.  (IV- 1)  for  those 
values  of  f  that  do  not  differ  by  more  than  approximately  0.001  from 
the  more  accurate  expression  for     f. 

By  assuming  that   b   (the  rate  of  boundary  movement)  has  been 
determined,    a  corrected  value  of  the  fractional  release  may  be  ob- 
tained by  subtracting  the  quantity    —    ( — J    from  the  actual  fractional 
release.     This  corrected  value  of   f   may  then  be  used  with  Eq.(III-22) 
to  determine    D. 

E.      Slab  Case  - -Moving  Boundary 

1        l/2 
For  large  values  of    — •  (3r    '     ,     the  fractional  release  is  approxi- 
mated by     (3r--the  percent  of  the  body  that  has  sublimed. 

STl/2       /ftTl/2  \ 

For  small  values  of    *- ,  (  *- <  0.2  )  ,     a  simplified  expres- 
sion for  the  fractional  release  may  be  developed  in  a  manner  similar 
to  that  for  the   spherical  case.      This  is  written  as 
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(  1)     P  =        0  (no  evaporation)        (4)     p  =     200 
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Fig.      IV.  1.       Simplified  expression  for  fractional  release 
for  the  sphere. 
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n       i  x 

2 


Hence,    D    may  be  determined  in  a  manner  similar  to  that  for  the 
spherical  case.      For  small  periods  of  time,    the  correction  factor 
y    Pt     or   —  ( —  J       is  negligible;  thus  a  straight-line  plot  may  be  ex- 
pected; and  over  longer  time  intervals  the  plot  may  be  expected  to 
show  an  increasing  slope. 

F.      Some  Numerical  Considerations 


Using  the  values  of  Shaked     for  his  sample   #  108  at  2040°C,    one 


obtains 


'6 


t'/z=  I  50  (sec)'*  (t=Z25*/0*sec); 


and  OL  -  0.2  5  CY*)> 

Therefore,         ?=  ^     -  2  3 '8  «l$*         (t* '  =  I \  5 7 '* 10* ') 


The  expected  fractional  release,    with  the  use  of  Eq.  (111-22),    would  be 

Shaked  estimated  that  the  value  of   b   he  encountered  during  the 
anneal  of  this   sample  was  0.3  u/h.      Because  Shaked's   sample  was  an- 
nealed in  a  Knudsen -type  cell,    this  evaporation  rate  was  considerably 
lower  than  would  have  occurred  if  the   sample   were  completely  exposed 
to   vacuum. 

Thus. 


(3-%^=  3/6, 


and,    if  the  value  of    f    calculated  above  is  substituted  into  Eq.  (IV- 1), 
one  has 
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with    D    representing  a  modified  value  of  the  diffusion  constant 

D  =  ¥.  I *  lo~ ,z  c«%€C  , 

which  is  61%  less  than  that  calculated  by  neglecting  evaporation. 

If  the   sample  had  been  completely  exposed  to  vacuum  at  this 
temperature,    the  evaporation  rate,    b,    estimated  from  thermodynamic 
data,      would  be  about  150  \i/h,    from  which 


j3  =  156,000- 


4 
For  an  anneal  time  of  2.25  X  10     sec  and  an  actual    D   of 

-  \  2  2 

6.6  y  10  cm  /sec,    it  follows  that 

'A 


Using  these  data  in  Eq.  (111-2  1),    one  predicts  a  fractional   release  of 
1-  0.757.      Now,     if    the  experiment  had  been  performed  in  this  way, 
and    if   the  effects  of  evaporation  were  ignored,    these  values  of    f,    t, 
and    a    used  in  the  nonevaporation  Eq.  (Ill -22)  would  yield  a  value  for 
D    of 


D-133  x  io~a  CYr*Z 


/sec, 

which  is  approximately  21,000  times  as  large  as  the  actual  diffusion 
coefficient  for  this  temperature. 

These  calculations  illustrate  the  importance  of  either  mimmiz 
ing  evaporation  experimentally,    or  of  taking  it  into  account  analyti- 
cally whenever  one  obtains  diffusion  data  at  high  temperatures. 
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V.     CONCLUSIONS 

An  equation  has  been  developed  that  relates  the  fractional  re- 
lease,   f,    of  a  tracer  element  in  a  geometrical  body  under  conditions 
of  evaporation  of  the  body  and  diffusion  of  the  tracer  through  the  mov- 
ing boundary  of  the  body,    with  these  conditions  occurring  simultane- 
ously.     This  may  be  better  expressed  as  the  fractional  release  of  the 
tracer  element  through  a  moving  boundary. 

For  small  values  of  the  quantity    (3,    the   release  curve  is  approxi- 
mately a  straight  line  whose   slope  is   slightly  greater  than  that  for  the 
stationary -boundary  case.     As    (3   becomes  larger  in  magnitude,    the 
slope  of  the  curve  increases  with  increase  in  time  and  approaches  a 
limiting  value    (see  Fig.  III. 2)    for  the  values  of    (3    and    t    considered. 

It  can  be  seen  that  to  completely  ignore  the  evaporation  of  the 

body  may  lead  to  an  incorrect  value  for  the  diffusion  constant.      The 

1/2 
curve  of  f   vs   (t)    '       can  appear  to  be  a  straight  line  when  the  bound- 
ary motion  is  not  zero,    but  the   slope  of  this  line  is  not  the  same  as 
that  for  the  zero -boundary-motion  line;  hence,    different  values  of   D 
would  be  obtained  unless  a  correction  were  applied  to  the  moving- 
boundary  case. 

Simplified  expressions  have  been  developed  that,    for  small 

l/2 
values  of  the  variable    (3t         /Z,     may  be  used  to  make  a  more  accurate 

determination  of  the  diffusion  constant,     D.      With  these  expressions 
one  can  predict  a  release  curve  whose  slope,    when    f    is  plotted 
against  the  one -half  power  of  time,    will  increase  with  time.      The  cor- 
rection made  to  the  fractional -release  equation  for  the   stationary - 
boundary  case  is  independent  of   D,    and  by  subtracting  this  correction 
from  the  actual  release,    one   should  obtain  a  line  of  constant  slope. 

It  also  should  be  possible  to  fit  experimentally  obtained 
fractional -release  curves  to  the  curves  obtained  in  Sections  II  and  III 
of  this  report,    and  thus  to  determine  both  the   rate  of  boundary  motion, 
b,    and  the  diffusion  constant,    D. 

Experimental  verification  of  the   simplified  and  the  approximate 
relations  for  the  fractional  release   remains  to  be  demonstrated. 
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APPENDIX 

A.      "Conventional"   Analysis  of  the  Slab  Problem 

Solutions  to  the  heat-conduction  equation  for  a  slab     -a<x<a, 
with  a  constant  initial  temperature    Tn,    and  with    T  =  0    at   x  =  ±a    are 
given  as 


-r/v+\-    ¥7Z  5  ('if  exp  r  ten+lfiftl        fa+Qirx 


(A-l) 


and 


TtiihT -TSC-flerfr  (P-^k-x  +etx.0n+l)a.+xl  . 
H*ty    i0    lo^t  V  1erTC-2(Jtt-)*.    +e^~J(E^j      (A"2) 

7 
The  average  temperature  in  the   slab  is  given  as 

T  -  SJL  Z  U—  exo  F-  £&2t£z?£  7  (A.3) 


and 


Using  an  analogy  to  the  heat-conduction  solution,    it  will  be   shown  that 
the  concentration  of  a  tracer  element    C(x,   t)    for  any  time     t    and  any 
position   x,    where     -a<x<a,      the  initial  condition    C(x,   0)  =  C    ,    where 
C_    is  a  constant,    is 


and 


CkthC.-ClW{**&ffi+erG%fa*y  . 


A-6) 
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Here,    D    is  the  diffusion  coefficient  for  the  particular  element  being 
studied.      The  corresponding  fractional  release,    f,    defined  as  the  ratio 
of  the  decrease  in  average  tracer  concentration  to  the  initial  average 
concentration,    is 


and 


■f-  zmfyA^i^^c  ftp] .    (A.8) 

A.l.      Tracer  Concentration  —  Exponential  Form 

By  using  the  method  of  separation  of  variables,    it  may  be 

o 

shown     that  a  solution  for  the  diffusion  equation 

and  (2    (±a,t)   =    0    y  (A'9b) 

may  be  expressed  as      C    (x,   t)  =  F   (x)  G  (t)  ,  (A-10) 

n  n  n  ' 

where  F  (x)  =  A     cos  k   x  +  B      sin  k   x  (A-lOa) 

n     '  n  n  n  n  ' 

and  G„(~t)~<^eXp(-k»l?tj.  (A-lOb) 


Applying  the  boundary  conditions    (A-9b),    yields    B     =0    and 

k     =—  (2n+l),    where    n=0,    1,  2,  •  •  •  .      Applying  the  initial  condition 

(A-9a)    yields    a=C_,    so  Eq.  (A-10)    becomes 
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C^t)  =  C0  K  cos  K  x  eyp  (-&X) 

The  coefficient   A      may  be  evaluated  by  setting   t  =  0   and  using  the 

n  9 

property  of  a  Fourier  series     that 

ft 


thus,    /\n=±fr(x)cos£e&Zdx   • 

-a 

This  may  be  integrated  directly  to  yield 


where    f(x)   =   1.       (A-13) 


4f-i)n 

A     =      l\     V  n  =  0,  1,   2  •  •  •   .  (A-14) 

n    IT  (2ml) 

Therefore,    the   substitution  of  Eq„  (A-14)  into  Eq.  (A-ll)  gives 


y  .    (A-i5) 


L"lW    rtz^T)  €X^L      ice-       r5  2a 


Any   sum  of  solutions  is  also  a  solution,    so  the  final   result  is 

my-  f*|45  *P[-&$^]cos&^  •  «A-i6, 


For  long  periods  of  time,    the  tracer  concentration  should  become 
zero  as     t  -*■  <x>  .      By  noting  that  the  exponent  in  each  term  approaches 
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-oo    as    t  -»•  oo,    it  is  clear  that  the  tracer  concentration  will  indeed  ap- 
proach zero  as  a  limit. 

A. 2.      Fractional  Release — Exponential  Form 

The  fractional  release,    f,    may  be  determined  from 

_f_     Qo-Qf^)  (A-17) 

where    Q(t)    is  the  total  number  of  tracer  atoms  per  unit  area  contained 
in  the  solid  at  any  time    t,    and   Qn  =  2aCn    is  the  total  number  of  tracer 
atoms  per  unit  area  contained  in  the   solid  at  time    t  =  0: 

r      2aC0  -fcfa)dx 

-a. 


For    t  -*■  oo,    the   summation  term  approaches    0    as  a  limit;  hence,    the 
fractional   release  approaches    1    as  a  limit.      For    t  =  0,    Eq.  (A-18)  be- 
come s 

-C  -     \—&-    5^  -, l—r-    ■  (A-19) 

oo 

111                         ^~ •             1  9 

The   sum    —  +  —  +  —  +  •  •  •  =     2_,     is   shown     to  be  equal  to 

l2       32       52  n  =  0   (2n+l)2 

TT2  8  TT2 

— — ,    so    f  =  1  -  —  X  —  =  i_i  =  o    at  time     t  =  0,     which  is  as  expected. 
8  2         8  ^ 

TT 
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Thus,    it  has  been  demonstrated  how  Eqs.  (A-5)  and  (A-7)  were 
obtained,    and  it  has  further  been  shown  that  the  tracer  concentration 
and  the  fractional  release  approach  the  proper  limits  for  large  periods 
of  time,    and  at  time     t  =  0. 

It  will  now  be   shown,    with  the  use  of  La  Place  transforms,    that 
Eqs.  (A-6)  and  (A-8)  are  the  solutions  of  Eqs.  (A-9)  and  (A-17),    re- 
spectively. 

A.3.      Tracer  Concentration  —  Error -Function  Form 


Let  ^Cfx.t)    _    J,    M 


(A-9) 


with       d  (xto)  =  C0  (A"9a) 

and  £  (±  aft)  =    0  J  ( A-9b) 

by  symmetry,  ^C(  OjC )  _     s\    t  (A-20) 

Now  the  La  Place  transform  of    C(x,  t)  =  C(x,  p); 

o 

Making  the  transformations  of  Eqs.  (A-9),    (A-9b),    and  (A-20),    one  ob- 
tained 

with  C(±atP)    —     0  (A-22a) 

and  ^C(Q,p)     _    q      .  (A_22b) 

dx 
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-> 
Defining    q    =  p/D    and  by  substituting  in  Eq.  (A-2Z),    one  obtains 

6>X^  ^  ^ 

Solving  Eq.  (A-24)  for     C(x,  p),     one  has 

C(x,p)  =  Aef-t~Be~f+  %■  ■  (A-"> 


Upon  applying  Eq.  (A-22a)  and  Eq.  (A-22b)  to  Eq.  (A-25),    it  follows 
that 

C(W-      p  (ef+eTj  +  p  (A"26) 

Multiplying  the  first  term  top  and  bottom  by    e    "    ,     one   sees  that 

C(KP) -  -^-[ef^e-t^e-^  j* ,      (A_27) 

and  upon  expanding    (1  +e      ^    )  by  the  binomial  theorem,    one  finds 

l+e-2f ]"'=  I-  e*T+(e-*fT-(e-*r? 


-4—      i  /   / 

(A-28) 


oo 


By  substituting  (A-28)  into  Eq.  (A-27),    it  follows  that 


<*-*>  ^-<k(*+afl  a.  Qi 


c(*lP)=  -%ji(-tfe-*T  [ene-H+  £ 


(A-29) 
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Taking  the  inverse  transforms,    one  obtains 

For  long  periods  of  time,    the  tracer  concentration  should  be- 
come zero  as     t  -»■  «,      If   t   is  allowed  to  approach  infinity  in  Eq.  (A-30), 
the     erfc  —  -*   1.0,    and  Eq.  (A-30)  becomes 


o=> 


I  in   C(x,t)=  Co-C0ZH)"(Z) 

(A-31) 

=  C0-C0(Z-2+2-2  +  -). 

which  appears  to  be  nonconvergent  and  equals    C.    or    -C_,    depending 
upon  whether  the  last  term  of  the  expansion  is  negative  or  positive. 
This  is  a  periodic  function,    however;  and  if  the  function   f(x)    is  defined 
as 

f(x)  =  2  if  O^x^tt,  (A-32a) 

and  f(x)  =  0  if  -tt  «  x  <  0  ,  (A-32b) 

then  the  coefficients  of  the   Fourier   series  for    f(x)    are 

—  TT 

Qn=-jfj2C0S»xdx   =£}   &»    W]J    =    0   >  (A-33b) 


and 


■7r  ,TT 


k-^rJZsm  m  dx  =  ^#Cos  nx]   =  ^(/-Cosmr)  ,  (a-33c) 


and  the   Fourier   series  for    f(x)    becomes 

_    ,      J£  f  *sin  x        sty)  5x    .   sin  5x      , . , ) 
~    I  +  rr  I       I  3  5  J 
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-     /  +  -L    ^ -^~r    Sir)  (Zn  I  l)X   •  (A-34) 

9 
On  page   180  of  Sokolonikoff  and  Redheffer,      it  may  be  seen  that    f  (x) 

is  well-defined  on  the  closed  interval     -tt  <  x  <  tt    except  for  the  simple 

discontinuity  at    x  =  0,      has  a  finite  number  of  maxima  and  minima, 

and  is  defined  for  other  values  of    x    by  the  periodicity  condition 

f  (x  +  2tt)  =  f  (x) ;    therefore,    the  Dirichlet  conditions  hold  for    f(x),      and 

the  Fourier  series  for     f(x)    converges  to 

-f(x)  =  -£[-f(x+) +{(*')]  (A"35> 

for  every  value  of  x. 

In  particular,    at    x  =  0, 

-f(0)  =  i  \f(0+)  +-ff0-)]-£(2+o)  =  I  ■  (a-36) 

Because  the  coefficient  of  the  second  term  in  Eq.  (A-31)  is   shown  to 
converge  to  the  value    1,    the  tracer  concentration  will,    in  fact,    ap- 
proach zero  as  a  limit  as     t  -*■  <» . 

A. 4.      Fractional  Release — Error -Function  Form 

Calculating  the  fractional  release  from  Eq.  (A-17),    one  obtains 


c     ZaC0~  fC(x,-h)dx 

'  ZCL  C  o 
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which,    with  the  use  of  Eq.  (C-ll),    becomes 


(A- 


By  letting    n  +  1  =  m,      the  last  two  terms  may  be  combined,    and  by 

-1/2 
noting,    from  Eq.  (C-7),    that    ierfc  0  =  tt     '     ,      one  finds  Eq.  (A-38)  to  be 

-?=  z($ffr++i%Wi''*&}-      (A"39) 

1/2 

For     t  ->  0,      the  terms  in  Eq.  (A-39)  involving     ierfc  k/t   /      become  zero; 

thus,   only     the     first  term      j—    ( )   '  remains.      For     t  =  0,      this 

tt1/2    ^a2/ 
term  and  also  the  fractional    release  are   zero. 


For  large     t,     as     t  -*■  00 ,      the  terms  involving  the 
na/tDt)1' 2->  ierfc   0 
tions  Eq.  (A-39)  becomes 


A     I  J  _  A     I 'o 

ierfc  na/(  Dt)    '      -*■  ierfc   0,      where     ierfc   0  =  tt  ;     with  these   substitu- 


By  using  the   same  analysis  as  that  for  the  tracer  concentration,    it  can 
be   shown  that  the   summation  term  converges  to    -1;    thus  Eq.  (A-40)  is 


-44 


which  is  indeterminate;  therefore,    it  is  necessary  to  refer  to  Eq.  (A-37) 
to  resolve  the  indeterminateness .     H^^CQ. 


(A-41) 


=  22(-0- 


This  previously  was  shown  to  be   1,    so 

llW     T~    I    '  (A-42) 

It  has  now  been  shown  how  Eqs.  (A-6)  and  (A-8)  were  obtained,    and  it 
has  further  been  demonstrated  that  the  tracer  concentration  and  the 
fractional  release  approach  the  proper  limits  for  large  periods  of 
time,    and  at  time     t  =  0. 

A. 5.      Correlation  Between  the  Two  Forms  for  Fractional  Release 

Using  the  method  of  Booth,      one  can  show  that  Eqs.  (A-7)  and 
(A-8)  are  identical  for  values  of     <   1.      Here, 


_2  r~)t 

for     <   1.      Equation  (A-8)  becomes 

a2 


,ii/2  W 
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Rewriting  Eq.  (A-7),    one  has 


The   sum  ^    ^~^  X     -    /  -/-  g  <£  £>-™  *  (A-44) 

OO  o 

/"XU 
e  du.      Thus, 


&G 


'e'^du.  =   (%/*•  <a-«» 


where   x   for  this  case  is     ,    K    is   1  or  4,    such  that 

4a2 

00  A 


+Z2trn*  ~-{VxT- 


Now,    if  both  sides  of  Eq.  (A-46)  are  integrated  over  the  limits  0  to  x, 
the  result  is 

m-i  "I  ' 

which  may  be  arranged 

Upon  substituting  Eq.  (A-48)  into  Eq.  (A-43),    one  finds, 


■46 


=  /  +  4 


77' 


G? 


¥C? 


</z 


£  IJL  fAdgtY^    £z±££    t_l 


77- ■ 


V 


4  a? 


3a3 


=  /+&(#)*- ££ 


77 'M  ce 


77"^  n=i 


n- 


n=l    n< 


f- 


f1  1  9  TT 

The        2_,     —     is   shown     as  equal  to    —  ;    hence, 


(A-49) 


(A-50) 


For  values  of 


7T2Dt 

a2 
.2 


<    1,      Eq.  (A-8)  reduces  to  Eq.  (A-50);  thus,    for 


these  values  of     ,      Eqs.  (A-7)  and  (A-8)  are  equivalent. 

a2 

Solutions  of  the  form  of  Eq.  (A-7)  converge  slowly  for  small 
-> 
values  of    Dt/a    ,     but  solutions  of  the  form  (A-8)  are  rapidly  conver- 

7 
gent  for  the   same   small  values.        This  statement  also  applies  to  the 

equations  for  the  tracer  concentration,    i.  e.  ,    to  Eqs.  (A-5)  and  (A-6). 
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B.      "Conventional"    Analysis  of  the  Sphere  Problem 

Solutions  to  the  heat-conduction  equation  for  a  sphere,     0  <  r  <  a, 

with  a  constant  initial  temperature    Tn,    with    T  =  0   at    r  =  a,    and  with    T 

7 
finite  at    r  =  0,    are  given  as 

T(0t)  =  -2T0g(-ir*xp(-*g£-)  «  <*--> 

n-i  / 

and  T(rt)=T0-^g[erJc^-er£%$$]  .    (.4) 
T(o,i)  -T.-$pl£*P  tWJ  ■ 


[See  note  following  Eq.  (B-6a).] 

7 
The  average  temperature  in  the   sphere  is  given  as 


lav  ~  ~W^  fy,7r*    Xr\       *x      / 


and 


(B-3) 


-&  =Z-6T(g)Y^M^$q+3ri(g).  (b-4, 

Using  an  analogy  to  the  heat -conduction  solution,    one  can  show  that 
the  concentration  of  a  tracer  element   C(r,  t)    for  any  time     t    and  any 
position     r,    where     0  <   r  <  a    and  the  initial  condition    C(r,   0)  =  C», 
with    Cn    a  constant  and  with    C(0,   t)    finite,    is 
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C^t)  =  -^^tfeKp(-Jf^2L)  „»  Ml  ,  (B-s, 


and 


(B-6) 


(B-6a) 


[it  will  be  noted  that  Eq.  (B-6a)  differs  from  ( B-Za)  - -other  than 
in  the  change  of  notation- -by  a  factor  of  two  in  the  summation  term. 
It  will  subsequently  be   shown  that  Eq.  (B-6a)  is  the  proper  expression 
for  the  tracer  concentration  at  the  center  of  the   sphere.] 

The  corresponding  equations  for  the  fractional  release  of  the 
tracer  element  are 


and 


B.l.      Tracer  Concentration  —  Exponential  Form 
The   radial -diffusion  equation 

^(rtt)  =  ^  pX  (n-t)      2   )C(r,t) 
at  L    dr*  r       dr 


(B-9) 


-49- 


may  be  reduced  to  the  linera -diffusion  equation  (A-9)  by  making  the 
substitution    u(r,  t)  =  rC(r,  t) 

with  the  following  boundary  and  initial  conditions 

u(r,   0)  =  rCQ  ,  (B-lOa) 

u(a,  t)  =  0         ,  (B-lOb) 

and  u(0,  t)  =  0         .  (B-lOc) 

A  solution  for  Eq.  (B-10)  may  be  obtained  by  using  the  method  of 
a 
separation  of  variables 

U   (r,  t)  =  F   (r)  G  (t)  ,  (B-ll) 

n  n  n 

where  F   (r)  =  A     cos  k    r  +  B     sin  k    r  ,  (B-lla) 

n  n  n  n  n 

-k  2  Dt 

and  G  (t)  =  a  e     n         .  (B-llb) 

n 

Applying  the  boundary  condition  (B-lOc)  yields    A     =0.      The  boundary 

condition  (B-lOb)  gives   k     =  ,     whereas  Eq.  (B-lOa)  gives     a  =Cn, 

(n=  1,   2,   3  •  •  *  )»    and  Eq.  (B-ll)  becomes 

U„(nt)  =  C0Bn sin  k„r  exp (-krzTJt) 

=  C0Br,*nWexp(-£g2£). 


The  coefficient    B       may  be  evaluated  by  setting    t  =  0    and  by  using  the 

n  9 

property  of  a  Fourier   series  that 

i    / 


■50 


Thus, 


Bn  =  -fi-ft(r)  sin  -Zf-  dr  >        wh< 


iere  f (r)  =  r  .  (B-14) 

-a. 


This  may  be  integrated  by  parts  to  give 

~R^-—(-))n         .  (n=  1,  2,  3,---)-  (B-15) 

The   substitution  of  Eq.  (B-15)  into  Eq.  (B-12)  results  in 

Ujrfl  =  -2^hTexp  f-t&gp£]  s*  J^rt  .        <B-i6) 

Any  sum  of  solutions  is  also  a  solution;   so  the  final  result,    upon  using 
C(  r,   t)  =  —  u  (  r,   t),      is 

and 

CfothPZ  C(nt)  =  -2C0ic-/fexp[-^^]  -       (B-m, 


where  the  relation 

lim     sin  x 
x-*0    ~x 


1 


has  been  used. 

For  long  periods  of  time,    the  tracer  concentration  should  be- 
come zero  as     t  -»■  <» .      By  noting  that  the  exponent  in  each  term  in  both 
Eqs.  (B-17)  and  (B-17a)  approaches     -oo     as     t  -+  °o,      it  is  clear  that 
the  tracer  concentration  will  indeed  approach  zero  as  a  limit. 
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B.2.      Fractional  Release — Exponential  Form 

The  fractional  release  may  be  determined  from 


■f 


_  Q0-Q(V  ,  (A.17J 


where    Q(t)    is  the  total  number  of  tracer  atoms  contained  in  the   solid 

43 
at  any  time     t,      and    Qn  =  —  Tra   Cn    is  the  total  number  of  tracer  atoms 


contained  in  the  solid  at  time     t  =  0.      Thus, 


f= 


=  i-M-^W^^Hv* 


(B-18) 


For     t  -»-oo,    the   summation  term  approaches  zero  as  a  limit;  hence, 
the  fractional  release  approaches   1  as  a  limit.      For    t  =  0,     Eq.  (B-18) 
becomes 


-f=  \-Ql£4s   ■  (B-19) 


2 


The   sum      — —  +  — —  +  — —  +  ■  •  •  =     Y*    is   shown     to  be  equal  to     ; 

l2        22        32  n=l  n2  6 

/  2 

so,        f  =   1 —  X  —  =1-1  =  0      at  time       t  =  0,      which  is  as  expected, 

2        6 

TT 

Thus,    it  has  been  demonstrated  how  Eq.  (B-5)  and  (B-7)  were 

obtained.      It  has  further  been  shown  that  the  tracer  concentration  and 

the  fractional  release  approach  the  proper  limits  for  long  periods  of 

time  and  at  time     t  =  0. 

7 
It  will  now  be   shown,    by  using  La  Place  transforms,      that 

Eqs.  (B-6)  and  (B-8)  are  the   solutions  of  Eqs.  (B-9)  and  (A-17),    re- 
spectively. 
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B.3.      Tracer  Concentration  —  Error-Function  Form 
We  have 


(B-9) 


with  C(a,   t)  =  0      ,  (B-9a) 

C(r,   0)  =  CQ  ,  (B-9b) 

and  v^S-J-J-Q  ,  by   symmetry  .  (B-20) 

The  La  Place  transform  of   C(r,  t)=C(r,  p)     is 


Making  the  transformations  of  Eqs.  (B-9),    (B-9a),    and  (B-20),    one  ob 
tains 


pC(rlF)-C0  =  V[l^^^], 


(B-22) 


where  C(a,  p)  =  0  ,  (B-22a) 

(B-22b) 


and  l^(0LP).  =    Q 

dr 


Defining 


q2   =  P/D  (B-23) 


and  substituting  Eq.  (B-23)  in  Eq.  (B-22),    one  has 


^^  +  2iC^.^crr,p)--§-    ■  (b-24) 


Solving  Eq.  (B-24)  for    C(r,  p),    one  obtains 


ecr,p>^er+f  e-r^. 
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Upon  applying  Eqs.  (B-2Za)  and  (B-22b)  to  Eq.  (B-25),    it  follows  that 

C(r(p)-      rp[e-y_e^j  +  P  <B-26> 


Then,    multiplying  the  first  term  of  Eq.  (B-26)  top  and  bottom  by 
e    4    ,     one  finds 

arlPy--^[e-ru^li-e-zr]-'+^,    (B.27) 

and  expanding    [l  -  e     q  by  the  binomial  theorem,    one  gets 


n-o 

Substituting  Eq.  (B-28)  into  Eq.  (B-27),    one   sees  that 

Upon  taking  the  inverse  transform,    one  finds 
and 

CCO^)=^oCrr(-t)  =  Co-#'  <B-30*> 

10 

which  is  indeterminate.      Applying  L'Hospital's   rule        to  the   second 

term  to  eliminate  the  indeterminateness,    one   sees  that 

x-^a  Q(x)        x-*a    g'  fx) 
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Now, 


Oji-r 


Cm-  C.-  >;Z  &  |0[er/c%^-  erfc%| 


^BUM-^^. 


As  an  additional  check  on  Eq.  (B-30b),    it  can  be   shown  that  the 
tracer  concentration  approaches  zero  as  a  limit  as     t  -*■  °o . 
The   sum 


ao 


/X3 


is,    for   small    x,    essentially  /     e  du, 

-  OO 


(A-44) 


oo 


where 


e^u  ■  fef 


(A-45) 


—  oo 


Equation  (B-30b)  may  be  written  as 


C  (othc^Me^-m^M-^ 


Combining  Eqs.    (A-44)  and  (A-45),    one  finds  that 


(B-31) 


c^> 


& 


n-i 


(B-32) 
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and  by  substituting  Eq.  (B-32)  into  Eq.  (B-31),    one  obtains 


lirrj 


I'm  2dCo 


z  \  a 


HTrl?t\lz 


L 


Z 


m,/zl 


=  Co  -  C0  =  0  • 

B.4.      Fractional -Release — Error -Function  Form 

Calculating  the  fractional  release  from  Eq.  (A- 17),    one  has 


(B-33) 


f- 


(B-34) 

r 


l-&(rC(nt)<lr 

=  -fsJo  *"&  Le r+c  20*)*     erfC  »  J  d'  - 

Now,  Eq.  (B-34)  may  be  integrated  by  parts,    by  using  Eq.  (C-ll)  and 
by  letting    u  =  r    and 


Z  C5€p7S 


which  yields 


1-  i  [«#  ^'^%fe  -  ^i^^t 


+agmf,erfc^^^(rt$^^M^  = 


~2Wty> 
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By  letting    n+  1  =  m,     the     ierfc     terms  may  be  combined;  the 

2  2 

i    erfc     terms  will  cancel  except  for  an    i    erfc     0     term.      Noting  that 

-  1  /?  7 

the     ierfc  0  = -n     '        and    i    erfc   0  =  l/4,    one  finds  that  Eq„  (B-35)  be- 
comes 


t^6(§)V^I^W-]-^^. 


(B-36) 


1/2 
For    t  -*■  0,      the  terms  in  Eq.  (B-36)  involving    ierfc  k/t   /        become 

zero.      For    t  =  0,    the   remaining  terms  are  zero;  thus,    the  fractional 

release  is  zero. 


If  the   substitution  of  Eq.  (C-7)  is  made  for     ierfc 
Eq.  (B-36),    the   resulting  expression  becomes 

f-6(m,/z 


na 


l/Dt 


m 


<£■'      L 


(B-37) 


For  large    t,    as     t 

into  Eq.  (B-37).      Hence 


oo ,     — -  -*■  0;    thus,    Eq.  (B-32)  may  be   substituted 


\\m  -£=  l'irr} 


t-9 


*-»« 


l/'ftl 


+  7T-"z[Jr(rrDtf-,])-3(j$)] 


—  I  in? 


na. 


n^\ 


(Dt) 


Vz 
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This  is  equivalent  to     oo -»■  oo,    which  is  indeterminate.      Referring  to 
Eq.  (B-34),    one  finds  for  large  values  of   t   that  the  difference  be- 
tween the  two  complementary  error  functions  may  be  approximated  by 


2  [(2ml)d+r     fe"+'k-r] 


(2n+ifof\      2r 


r_fenW)V] 


which,    upon  substitution  in  Eq.  (B-34),    gives 

U  Jl  -  Urn    3  far*    f  e       L&2±!%£ldr 


-    1 1)77 


6 


t->«*  (£(ttM)^)q 


>WM-m-lM-W^ 


Applying  Eq,  (B-32),    one  finds 


I  in  S  -  li 

t  -=h>o  I         t- 


.a 


^f^{r%(^f-^^t]dr 


§5$kj}*  \_^^)k~  C^t)'/2  j  dr 


=  |  (>jr  =  #*3  =  /  • 


a3V  "'        3az 


Thus, 


lim 

t—OO 


f=  1 


(B-39) 


It  has  now  been  shown  how  Eqs,  (B-6)  and  (B-8)  were  obtained,    and  it 
has  further  been  demonstrated  that  the  tracer  concentration  and  the 
fractional   release  approach  the  proper  limits  for  large  periods  of 
time  and  at  time     t  =  0. 


Booth  has   shown  that  for  values  of 
(B-8)  are  identical. 


7T2Dt 


<    1,    Eqs.  (B-7)  and 
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Solutions  of  the  form  (B-7)  converge   slowly  for  small  values  of 

Dt/a    ,    but  solutions  of  the  form  (B-8)  are  rapidly  convergent  for  the 

7 
same   small  values.        This  statement  also  applies  to  the  equations  for 

the  tracer  concentrations,    i.  e.  ,    to  Eqs.  (B-5)  and  (B-6). 
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C.      Properties  of  the  Error  Function  and  Related  Functions 

7 
The  following  definitions  are  made: 

2        x      _z2 

erf  x  =   e  dz  ,  (C-l) 

\FJT  0 
erf   oo=   1,  (C-2) 

erf(-x)  =   -erf  x  ,  (C-3) 

2    r  -z2 

erfc    x  =  1   -  erf  x  =   /      e  dz  ,  (C-4) 

n/"tT     X 

ierfc   x  =     /      erfc    z   dz  ,  (C-5) 

x 


00 

i    erfc    x  =     /      i  erfc    z    dz  ,  (C-6) 

x 

where  i    erfc    x  =  erfc    x  ,  (C-6a) 

2 
ierfc    x  =    -±—  e  -  x  X  erfc    x  ,  (C-7) 

nTtt 


7  \ 

and  i   erfc    x  =   —   [erfc    x  -  2x   ierfc   x]    .  (C-8) 

4 
The  general  recursion  formula  is 

2n   inerfc    x  =  in~    erfc    x  -  2x    i"      erfc    x  .  (C-9) 

If  both  sides  of  Eq.  (C-6)  are  differentiated, 


d      .n+1      c  d        f   .n       .  , 

—    i         erfc    z  =    — -      /     i    erfc    y    dy 
dx  dx    J 


.n      ,  dz 

=   -l    erfc    z    . 

dx 


(C-10) 
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Integrating  both  sides  o£(C-10)  and  rearranging  yields 


f.n 
J   1    er 


fc    z    dx  = 


1       .n+1 

i  erfc    z 


dz 
dx 


(C-ll) 


dz 
provided    —    is  a  constant  and  the  constant  of  integration  is  neglected; 

dx  ....  n 

then  (C-ll)  is  the  indefinite  integral  of    i    erfc   z  dx,     where    z    is  a 

function  of   x. 

The  values  of   i    erfc(-x)    will  now  be   studied  to  obtain  an  inver- 
sion formula,    so  that  the  tabulated  values  of  the  positive  argument 
may  be  used. 

erf ( -x)  =   -erf  x 


Noting 
and 
one  obtains 


erfc    x  =   1   -  erf   x  , 

erfc  (  -x)   =    1   -  erf  (  -x)   =    1  +  erf   x 

=   1  +  (  1  -  erfc    x) 

-  Z   -  erfc    x  . 


(C-3) 

(C-4) 

(C-12) 


Thus, 


,                   1        -x  , 

lerfc    x  =   e  -  x   erfc    x 


nTtt" 


(C-7) 


and 


ierfc(-x)   =   e 

N  7T 


(  -x)erfc  (  -x) 


e  +  xerfc  (  -x) 


\Ttt- 


n/tt 
=  2x  + 


-  e  +  x(  2  -  erfc    x) 


1        -x' 
—  e 


xerfc    x 


■] 


=  2x4-  ierfc    x 


(C-13) 


So, 


i    erfc    x  =    —    [erfc    x   -  2x    ierfc    x]  , 

4 


(C-8) 


(C-14) 
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2  1 

and  i    erfc  (  -x)  =   —   [erfc  (  -x)   -2  (  -x)  ierfc  (  -x)] 

4 

-   —   [  ( 2  -  e  rfcx)  +  2x  (  2x  +  ie  rfc    x)  ] 
4 

=   —  -   —  (erfc    x  -  2x   ierfc    x)  +  x 
2         4 

1    ,      2       .2       , 
=   —  +  x      -  i   erfc    x  . 

2 

By  using  the  general  recursion  formula  (C-9)  and  the  relations  just 
developed,    it  may  be  shown  that  succeeding  terms  are 

3  113         3 

i    erfc  ( -x)  =   —  x  +   —  x     +  i   erfc   x  ,  (C-15) 

2  3 

i   erfc  ( -x)  =   —  +  —  x     +   —  x      -  i   erfc    x  ,  (C-16) 

16         4  12 

and  i   e  rfc  (  -x)  =   — -  x  +   —  x     +   —  x     +  i   e  rfc    x  .  ( C  -  17) 

16  12  60 

From  the  foregoing,    a  general  inversion  formula  for    i   erfc(-x)    may 

be  obtained: 

.  .                             _    n  n-2 

.n       .     ,       .        ,      ,,n+l  .n       -           .     2x       ,       x 
l    erfc  (  -x)  =  (  -1)  i    erfc  x  +   + 


n!         2(n-2)! 

(C-18) 
n-4  _    n-m 


l6(n-4)!  (4-8  -12---  2m)(n-m)! 


where  n  -  m  =  1  for  n  =  3,   5,   7,  •  •  • , 

and  n   -  m  =  0  for  n  =  2,4,6,'#-. 
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D.      Solution  of  the  Slab  Problem  for  a  Moving  Boundary 
D.l.      Concentration  of  the  Tracer 

An  expression  for  the  fractional  release  of  a  tracer  element 

,  5 
through  the  moving  boundary  of  a    slab  will  be  obtained.      Chambre 

has  derived  a  general  solution  for  the  tracer  concentration  for  a  slab 

-(a  -  bt)  <  x  <  a  -  bt    and    0  <  t  <  — ,      where     b    is  the  rate  of  movement 

b 

,0. 


of  the  boundary 

This  equation  is  the   solution  of  the  diffusion  equation 


(D-l) 


(D-la) 


(D-2) 


subject  to  the   restricting  conditions 

C(X  ,<>)-{(*)     for 

K(0£)    -     fS  for 

and  C(«-^t,t)=     0  f°r 


0  <x  <&, 

(D-2a) 

0<t<t, 

(D-2b) 

(D-2c) 

For  this  problem,    the  initial  condition  (D-2a)   shall  be  assumed  to  be 


a  constant,    C.-.;  thus, 


C&o)  =  Cd      for    0<,x,<^ 


(D-3) 


If  this  substitution  and   S(x,  t;    £  +  2na)    are  made  into  Eq.  (D-l),    one 


has 


~h 
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~(K 


m-^/W"6*^ 


.     (D-4) 


This  may  be  simplified  somewhat  by  defining 

y  =  |  +  2na 


and 


z  = 


upon  which  Eq.  (D-4)  becomes 


CM)  = 


_  c 


2  \T5t 


(D-5) 


(D-6) 


2(TTDt)A 


-2(ptf    exp(-z*)d-. 


M^ 


(D-7) 


a.(2rt+l)  -. 


#(2n+/) 
■f 

The  first  term  may  now  be  integrated  directly,    and  if  exponentials  are 
combined,    one  has 


2< 


(D-8) 


The  exponentials  may  be   simplified, 


M-^=wt^rf~Hn^ 


=  7^{f-2tf(x+2")'t)+(x+?"btT+^-(*+2»6tf 
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_     / 


V2>t 


fy+(*-2nl>t)]Z+^(x-nbt)  ; 


W)-M*^-^*S 


2(Dt)J 


M^^HM1*^/^ 


(D-9a) 


(D-9b) 


^S-'W-#'"»*.|/«p(-ir^21}^  ■ 


&&»-l) 


Simple   substitutions  similar  to  Eq.  (D-6)  allow  the  last  two  terms  of 
Eq.  (D-10)  to  be  integrated: 
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(D-ll) 


_l_  v-Q 


Using  Eq.  (C-3)  and  rearranging  terms,    one  has 

C(x,t)=  &-\erf 

Co 


a+x    ,  a  r  a-x 


4^|r-,rexP[-^K-^x)][er-r^|-^;rx 

en — zvm^— 

With  the  use  of  Eq.  (C  -4),    i.e.,    erfc  x  =  1  -  erf  x,      one  finds 

C(x,t)  =  Co-  Mer/c^+er^*w* 


(D-12) 


+ 


■J-|(-/)"exP[-#^-^-#^  ^g^ 


ertc    %7^v^ — 


Dt)' 


(D-13) 


^£(->rexp[-^-^4[erfc  ^jfegfr"* 


-f-  Co 

z 


For    t  =  0,    because  erfc    oo  =  0,    Eq.  (D-13)  becomes 


C(x,   0)   =  C 


0 
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For    t  =  a/b,     which  corresponds  to  total  evaporation,    the  tracer  con- 
centration should  become  zero.      Hence, 


a^c-Mer^w^er/c 


a-x 


(D-14) 


+  T^K-^t^-f^-  erftfg 


f 


Since     |  x  |  <  a  -  bt    at    t  =  a/b,     x   in  Eq.  (D-14)  must  be  zero.      Thus, 

-=  Cb-Urfc[U¥f\  >  C0^!)f*mrj-crfc[i(ifl 


with  the  use  of  Eq„  (C-12), 


=Co[i-eAi(&]+2e&ir[i-rfc4faf 


(D-15) 


It  was  previously  shown  in  Eq,  (A-41)  that    2     )_,    ("1)      =  1»    similarly, 


n  =  0 


it  can  be   shown  that    2     V    (  -1)     =  -1.      With  this  substitution, 

n=l 
Eq.  (D-15)  is 


C(*frCt!-er&i(%fl-C0\j-erftifej'*]  =0-    <°-">     I 
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If,    in  Eq.  (D-13),     b  =  0,     which  corresponds  to  no  evaporation,    then 


(D-17) 


If    2n-l  =  2m+l,     Eq„  (  D- 17)  beco 


mes 


5" 


=  n  __  Co 

L°       2 


5V_ \f  fr^(2n+J)0L+X    ,  ^jZ.(2n+lh-)C\ 


(D-18) 


-■&-£ 
2  ^ 


which  may  be   recognized  as  Eq.  (A-6). 

D„2.      Fractional  Release  of  the   Tracer 

The  fractional  release  may  be  computed  from  Eq,  (A- 17),    as 
follows : 


h 


(2„ 


(OL-kb) 

2a.0a-J  CM)dx 

-fa-U.) 

2a.C0 
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r(a-kt) 


^£6 


using  the   symmetry  relation,    one  obtains 

,  fa-it) 


^\-^rfC(x^dx-  (D-i9) 


With  the   substitution  of  Eq    (D-13)  for     C(x,  t)s     it  follows  that 


_  Prjfl  2r,Ca-kt)+a-hX  Irlv 

cnc     zfw^ — J 

With  the  use  of  Eq,  (C-ll)  in  the   second  integral,    Eq,  (D-20) 
become  s 

,.ft-(.t)      /-,t^        J 
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=f#/M 


fat 


2<x-kt 


rra-it> 


(D-21) 


-^^^y^(-^^^-er0^^- 


The   remaining  integrals  in  the  expression  (D-21)  for     f     may  be 
evaluated  by  integrating  by  parts,    as  follows.    Let 


U-exp(-^)  ana     dv--erfc2"f(^-ydx  ■      *»» 


du^xP(f)dx  and     y=  2(Dtfie&  2f{^tX 


The  result  may  be  written  as 


(D-22) 


To  integrate  the   second  integral  by  parts,    let 


u=exp(- 


and     dy^er^2nt^ta^  dx 


D 


2(Dt)v* 


then 


du=-gexpf$)<ix     -a     v=-2(Ktcer>fcZn%l$*£h* 
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The   result  is 


9 


(D-23) 


Upon  applying  the  limits  to  Eqs.  (D-22)  and  (D-23)  and  inserting  the  re- 
sults into  Eq.  (D-21),    one  obtains 

I  =  fat +(bt f  [■    r  jg_      ■    r  za-ht  1 
'     a  +^;  LLerf  2^2  ~ieHcHJ[5^- 


a? 


f(-lfexp[-f(ny^l%-,^ 


-i^^SggR  -exp[^)]^r£%|$^  ,d-24, 

+  lertc   20*)''*    ^TtT^F^  2(Dt)*  ^  ^CDty^J 

After  eliminating  equal  and  opposite  terms  and  combining  exponen- 
tials,   it  follows  that 
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(D-25) 


-t 


-t 


2(DtJ 


^fl^mX^T^pl-^M>^ey^^^ 


At    t  -  0,     each  term  of  Eq„  (D-25)  becomes  zero;  thus,    the  fractional 
release  is  zero.      At  time,      t  =  a/b,     which  corresponds  to  complete 
evaporation  of  the  solid,    the  fractional  release  should  be   1,    as  shown: 


-f- 


+ 


7TE 


'/zT. 


ler-fc 


a 


OclVz 


zffi 


Cerfc-if$* 


+ 


'D 


A 


\db  J  ntm=i 


^rm^t'r^my* 


(D-26) 


-l^er-fc 


a 


z(%y* +  &erfc-  2TW1 "  l"erfc  zm*- } 


=  1 


If,    in  Eq.  (D-25),     b     is   set  equal  to  zero,    which  corresponds  to  no 
evaporation,    then 


f= 


(i&n.^c  o-t^fr]+&f &*"[**  sit 


\  o^  ] 


cen 


(D-27) 


-LerfcMn;  +ierfc$A  -  Cer-fc^ 


a 


"(DtF     u " "  ^'TW7*      ~"  "  (oty*- 
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by  letting     n  -  1  =  p    and    n  +  1  =  m    in  Eq.  ( D-27),    one  obtains 

-1/2 
Noting  that    ierfc   0  =  tt  ,     and  combining  the     m    and    p     terms,    one 

finds 

which  may  be   recognized  as  Eq„  (A-8). 
D„3.      Conversion  to  Dimensionless  Form 


If  one  defines 


t  =  —  (D-30) 

.2 


ha 

and  (3  =    —  ,  (D-31) 

sothat  ^-^    "f"     ^  0<^O     ,  (D-32, 

where     Pt  =  0    corresponds  to     t  =  0,     and    (3t  =  1     corresponds  to  com- 
plete evaporation,    then  Eq.  (D-13)  becomes 


s 


-£(-l)nexP  i-nfifti-fiQ-k]  \  jer-fc- 


0-B*d')-$ 


;    i   x 


O   r~~l/; 


(D-33) 


/<0  n~t 


2.  Z  ''* 


Equation  (D-25)  becomes 

-f --  61  +  T'/z(ierfc-§^  -  Cerfc  j0k) 


c*2  fji'W 


\ / 


1       '        1  <-    1  1    -J-  z,z//z- 


n,w: 


(D-34) 


r„K--,  ,  ,    -    1  r    -'-  <  (.  "■ 


and  Eq.  (D-29)  becomes 


(D-Z5) 


D.4.      Approximate  Solution  for  the   Fractional  Release 

If  Eq.  (D-34)  is  expanded  to  the  first  few  terms  in    n,      then 


m~i 
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-  Z'^-O^zfe^O^erA  fg@  -  ^^  ^ 


# 


-  Z*&ir(6pzf"exp[-6?(l-lty(ter(i  4^f?  -  cWc  ^|« 


-  Z^LMe^p^-rf]^  4^  -  L»erfc  *£&. 


7T" 


+ 


Now,    if  the  restrictions    —  <  0,01     and    Bt  <  0„5     are  made,    then 

a^ 
Eq„  (D-36)  may  be  closely  approximated  by 


Inversion  formulae  for     i     erfc  ( -x)    have  been  developed  in  Appendix  C, 

If  these  are  applied  to  Eq„  (D-37)  and  expanded  in  terms  of   m,     it  can 

'3T1/2nx 


be   seen,    if  sufficient  terms  are  taken,    that  terms  of    I  — cancel 

m             Bt1/2  V      2       / 

and  only  the     i     erfc    *- terms  are  left.      Thereupon  the  approxi- 
mate expression  for  the  fractional  release,     f,    becomes 


$--^2Z\crfc^\lm^J^erfc^  •         <°-38> 
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E.      Solution  of  the  Sphere  Problem  for  a  Moving  Boundary 
E.l.      Concentration  of  the  Tracer 

In  this  section  an  expression  for  the  fractional  release  of  a 

tracer  element  through  the  moving  boundary  of  a  sphere  is  obtained. 

'  5 
Chambre      has  derived  a  general  solution  for  the  tracer  concentration 

for  a  sphere     0  <   r  <  a  -  bt    and     0  <  t  <  a/b,     where     b     is  the   rate  of 

movement  of  the  boundary 

a. 

r/\.    r  (r-r^l , 


^ fj%)&r[nDk(y^S(r,t}yZ-,^  -       (E-u 


/?=/! 


tf. 


W  STr^j)  =^[4#J  -  e*F[-%£]  ■ 


This  equation  is  the   solution  to  the  diffusion  equation 


(E-2) 


subject  to  the   restricting  conditions 

C(r,   0)   =  f(r)            for  0<r<a,  (E-2a) 

C(0,  t)        is  finite  ,  (E-2b) 

and                            C(a-bt,t)  =  0           for  0  <  t  <  a/b  .                          (E-2c) 

For  this  problem  the  initial  condition,  (E-2a),    shall  be  assumed  to  be 
a  constant,    C„;  here, 

C(r,   0)  =  CQ           for  0<r<a.  (E-3) 
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If  this  substitution  and    S(r,   t;  £  +  Zna)     are  made  into  Eq.  (E-l),  then 

.0. 


-« P{-L^^}}% 


The  first  term  may  be  integrated  by  parts,    by  letting 

^5 


a 


nd 


<*v=exp[-^h> 


if  the   second  term  is  rewritten,    then 


C(nt)  -  %&  £*§§*+  erfc^J-  fe/^C^  d$ 


The  exponentials  may  be   simplified,    as  follows 


kz^}2  ^  "ik  [r-2K^ycr24-^t5l 


-  jJL  W% tnos-  tfnktj  -2rr -  Vvar-h ¥/?2<?-hrzl 
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^\^2*(r-2n*.+2ryfy  -tnaM  f/fif+r*] 


ixt 


=  zn^xr  r^-2wa  ^n£tf -  2  r  (r-2na+2»bt)  +f 

iiyt  U  J  J  (E-6a) 


+  *tnar-lfn*aJL-r*:] 


and 


HDt 


ivt 


[jZ+2$(r+2n<i-2r>ht)  i-(r+ZMi-2r>bt)' 


(E-6b) 


-  (r-t2n6L~2h)5tf+  Hnar  +V/)V*r*J 

-  ynar-¥nza*-rzl 


I 


H!Jb 


[5+(r+2ho--2nhtyf-J3A[-r-2ria.+hbtJ. 
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Integrating  the   second  term  of  Eq.  (  E-5)  with  the  aid  of  Eq.  (C-ll), 
and  also  substituting  the  results  of  Eqs,  (E-6a)    and  (  E-6b),    one  finds 


CM^&Lerfctfo+ertc^] 


(E-7) 


'  -a. 


-     Co 


-a. 


The   remaining  integrals  may  be   solved  as  before  by  integrating  by 
parts.      Let  ,  .        ^ 

u~  5 


and  let       J v  =  exp  [-  L^figs&^f]  d 


c 


J 


and         OV-  €XP}—- 7-r+i-  ^       I  $T  '  successively;  then 

'   L  HTtt  J    J 


C(r,t)=  ^[et4cj^+erka+r 


2(Dt)** 


4 
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a. 


-*-  (E-8) 

a* 


-j^expp£fa»kt+$fer6  r+z%Z$£+5  %• 


Upon  integration  of  the  last  two  terms  of  (E-8),    one  has 

"''  .  „„r.  r+2na.-2nkt-a.l 

Wa/  ■>._/;  r-2na_+2nkt+cJ 

-  tenc    Z(Dtyk       J 


Lertc     2(Dbyh — J 


By  using  Eq„  (C-12),    one  obtains 


r   r-2no.+2nbt-CL       „   r   r-2n<x+2nht+CL 

eric    zttyty*       ^erfc     2(Dty''*- 
-4-er+c  Z(Dty^     +  /Z-ertc     z Cot)''*- 


(E-lOa) 


and  e)PfC   2p%  "    2-  er~&2(DE)'/*'  (E-lOb) 
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Using  Eq„  (C-13),    one  has 


Ler^  sk^  2[i&y<**0& 


■h£     a-r  (E-llb) 


The   substitution  of  Eqs.  (E-10)  and  (E-ll)  give 


■■£exp[-f(»o-~ntt+r)][erCc  Zn%^^r 


2f    n. 
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-  C  (Dt^Jg  exp  Fnh  (narnbt-f-j\\ierfc  2ng_-2nloi:+a-n 

Upon  cancelation  and  rearrangement,    one  has 


(E-12) 


C  (r,±)  =C0-^  [eric %&-  erfc  ^ 
-j^iexP[-nh(narnht-rj][efc  Ss^+^r 

+  /*r-Cr  2nfa.-tid-a.-l-l 

^ertc      Z(6t)<h      J 

,  atcr.  Zn  fa-tit) -Q-+r  1 

-t-erk    2(--<y; — j 

_cer/c^g-tt)-a-rj 


(E-13) 


a+r 


0' 


For     t  =  0,     since  the     erfc  oo   =  0,     Eq.  (E-13)  becomes     C(r,   0)   =  C 
For     t  =  a/b,    which  corresponds  to  total  evaporation,    the  tracer  con- 
centration becomes   zero.      (The  proof  will  not  be  given  here.) 

If,    in  Eq.  (E-13),    b  =  0,    which  corresponds  to  no  evaporation, 


then 


tfr,t)=C0-^[er£^-er£ 


OL+r 


•2(5$' 


Tx 
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'tit)''* 


If    2n-l  =  2m+l,      Eq.  (E-14)  becomes 

_  C0a-fw!/zripUr  a-r       -p^-   g-+r    -j 


(E-14) 


+  r{W)&0  lteric  Z(£y*     ieric  2(Dty* 
Zr   fc0  lenc  zittfa       €ffc   2(Dt)'k  J 

r  \  aV  &  Lterfc2f&fcy^  - cer^  2(mM 
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5  r.,r.M^r     ^^fen+Qz+n 


which  may  be   recognized  as  Eq.  (B-6). 

E.2.      Fractional  Release  of  the  Tracer 

The  fractional  release  may  be  computed  from  Eq.  (A- 17), 

O0-Q(t) 

f  =   — ,  (A-17) 

Q0 

where     O(t)     is  the  total  number  of  tracer  atoms  contained  in  the  solid 

4  3 

at  any  time     t,      and    0n  =  —    Tra    Cn    is  the  total  number  of  tracer  atoms 

contained  in  the   solid  at  time     t  =  0. 

(ark-t) 


kVith  the   substitution  of  Eq.  (E-13)  into  Eq.  (E-16),    one  finds 
(a-yt)  /a.y.) 

o  o 

^er^2^a-fet)-a-r]dr 
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<M>  .  -(.eric    z(bt)V,     jar 

__3_  f rY^Y^r^pP#Cm-^t+r)7f ce rfe  *" fer^a-H- 

6  .  r        '7>)n{l 


■f- 


The  first  integral  in  Eq.  (E-17)  may  be  handled  directly,    whereas  the 
second  and  last  integrals  may  be  handled  by  integrating  by  parts  twice, 
as  follows  with  the  use  of  Eq  ,  t  C  - 11) 

UL=  r  >  dv,  =  ^erfc7&jxfc  dr>     where   n* 0J /, 


du  =  dr,  and     K  =  2fDt)//fe6nerfc^^^» 
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or  dv^Cerfc^dr,  and  V^-ZiOtf'Ckr-fc^^  ■ 

Hence,        f  =    j  -  £  ^b^^Z(^]\(^t)ier{c^ 

+£{$ J* &&%**)  ^  fa**  (*-*)  i'erk  fgg 

/a-fct)  y    " 

+errr  2r,(G-ht)-a.+r  ~\a 
^3^yr|exp[-^^-^t-0][ce^c^ra-bt^a-r 

o 
0 

_  ferr   2hfa.-bt)-g+rlJr 


=  sm-  ?(£T+m* 
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fa-it)  2(  Dt)  *         J 

.  <w#  ccnc2j5^77^Jar 

The  remaining  integrals  in  the  expression  (E-18)  for  the  fractional  re- 
lease may  be  evaluated  by  integrating  by  parts,    as  follows: 

ix~  rexp5^C"a-ntt-i-)]. 

and     v-2(btfierfcz^-^0^  ■ 

Several  repetitions  may  be  required  before  the  pattern  that  develops 
becomes  obvious.      The  result  is 
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2(Dt)Vi^ 


2(Dt)' 


The   similar  integral,    which  would  have     dv  =  ierfc    — -^— — dr  , 

2  n/  Dt 
differs  from  Eq.  (E-19)  only  in  the  power  of  the     ierfc  x    term.      If    m 

were  replaced  by     m+  1,     the  result    [Eq,  (E-19)]    could  be  used  as  it  is 

The   remaining  form  of  the  integral  to  be  evaluated  by  parts  is 

(X  =  jr  expl-^Cnar-nbti-r)], 

and  ^-       ^(^v     .^m £_(b^>/z 


After  several  repetitions  the  result  may  be  recognized  as 


(E-20) 


Upon  applying  limits  to  Eqs.  (E-19)  and  (E-20),    the   resulting  terms  are 


inserted  into  Eq.  (E-18)  to  give 


a. 


+ 


■f 


+ 


+te&&£j$$-°-  ] 


yy>=z 


(E-21) 


t-terA&ZlfegbL 


tm- 


+ 


m=3 


L  ertc      2(DC)'^      J 
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as  the  solution  in  full  form  for  the  fractional  release  of  the  tracer 
element  from  a  spherical  body, 

At    t  =  0,    each  term  of  Eq,  (E-21)  becomes  zero;  thus  the  frac- 
tional release  is  zero.     At  time     t  =  a/b,    which  corresponds  to  com- 
plete evaporation  of  the  solid,    the  fractional  release  should  be   1,    as 
follows : 

{=3-3  +  1-  6(M*er&zt®'*- ^er/c  W*  ] 


m^z 


E-22' 


=  3-3+1  =/■ 

If,    in  Eq.  (E-21),     b    is   set  equal  to  zero,    which  corresponds  to  no 
evaporation,    then 


-f 
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=  3(&f[lerfii  0+terfc^l  (E.23) 

'  ^1  av^f  l^  '-  W~   '  -^r'  -  (frc/^  ^^c  (Dty/2  j 
'  Ma*>  ^LCe^7Dt^^^6e^cf^^^^rc-^gJ. 

By  letting     n  +  1  =  m    and    n  -  1  -  p,     one  obtains 


-91- 


-  ^DtY  To  ;nr£n  n-u*?  fpyC.  rna.  ,  0  <  ,vrZ.  Mr.   + <<  fPrJ^  _Pa 


(E-24; 


tSi3crfz-^-/z-i-l5er^c  0-oWcO] 


Because     ierfc   0  =  tt~    /        and     i    erfc   0=l/4,    Eq„  (E-24)  becomes 


r-Ci/Dt^r-rr-'k.,  n^fphfi.JjAr,l      ^[Ot\    ,  (E-Z5) 


which  may  be   recognized  as  Eq.  (B-8). 
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E.3.      Conversion  to  Dimensionless    Form 


~ifa3 ,exp  t-^"(,^~*lf  [erfc  Zh(!^>/i'~r/a- 


If  the  definitions  (D-30)  and  (D-31)  are  now  utilized,    along  with 
Eq,  (  D-32),    then  Eq.(E-13)  becomes 

wc^-ji-'-r^] 

_  Co  r'/z[:rrr.  \-rlo-     jprfc  i+r&  1  <E-26) 

-Cer-fc  ^jjgfeiz'h] 

Equation  (E-2  1)  becomes 

-l2Z*Ht3erfc^-tWcj^L)-  GZ(c2er&^ 
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+3r'/2(/-^)J:/j'w'r2«/5r^mexpr-^r«wr/-/3f)][<;"'erfc^^|s^ 

<  77  J        (E-27) 

and  Eq.  (E-25)  becomes 


-f =  G  T'/2[^'/x+2^  Cer^c^-SZ  •  <E"28> 


E.4.      Approximate  Solution  for  the   Fractional  Release 

Now,    if  the   restrictions     —  <  0.01     ar 

a^ 
Eq.  (E-27)  may  be  closely  approximated  by 


Now,    if  the   restrictions     —  <0.0l     and     3t<0.5     are  made,    then 

a^ 


f =  ^pz-^zf^zf+5Z'/z(i-f3z)  ierfc-^  - 


•94 


(E-29) 


m=2 


"7"  2 


-i2r%j6/r^-^^TeVr-^^-/r'-^c^rfc^#^, 


^=5 


of  which  the     n=  1    term  is  most  significant. 


Thus,   f  =  3p^-3(pzT+(ptfi-3T'^l-/*?)ler{c££ 


7a 


-fir^f%6rr/-^)^/2-/2r%L^^ 


kw+L,  vw-l  >  ar'h^'1 


^ZH^%^rW\^f  ^erfc  ^ 


+6z^rtm^f^^g^e^  ^ 


2 


'4a 


(E-30) 


-6z(\^i(-T(^\M^r\^rk.=^ 


M=2 


'Z 


-12  ?:3/2£c-/)wrm-2X4)m-t4^T"3^^  ^ 

/r?=3  <*-  <*~ 


'/* 
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Application  of  the  inversion  formulae  for     i     erfc  ( -x)  - -developed  in 

Appendix  C--to  Eq.  (E-30),    upon  expansion  to  a  sufficient  number  of 

(3Tl/2\x 


terms,    reveals  that  terms  in    (  -t- — j       cancel  each  other.      There- 
upon the  approximate  expression  for  the  fractional   release  will  reduce 
to 

f  =  3prc-3(pzT+((3i:?-3Z(i-pt)--3Z 
+  6  z\i-fiz)Lerfc  -g^V  12  Z(l-/ez)Lzerfc  $& 


+  5ZHl-r)i(«T(^7L^r-fc  f^ 


(E-31) 


+  eZ%(mHXVm(4fTlmerfc  -^ 

In  general,    the  last  three  summations  of  Eq.  (E-31)  may  be  neglected 
in  comparison  with  the   remaining  terms 

f  =  3fz(l-(5z)+(f5Zf  -  3Z-  3T(i-f3Z:) 
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NOMENCLATURE 

a  Characteristic  dimension  at  initial  conditions;  half  thickness  for 

a  slab;  radius  for  a  sphere  (cm) 

b  Rate  of  boundary  movement  (cm/sec) 

C  Concentration  of  the  tracer  element  (atoms/cm    ) 

Cn  Initial  concentration  of  the  tracer  (atoms/cm   ) 

D  Diffusion  coefficient  for  the  tracer  atoms  (cm   /sec) 

f  Fractional  release  of  the  tracer  element 

Q  Number  of  tracer  atoms 

QO  Initial  value  of  Q 

r  Space  coordinate  in  spherical  geometry  (cm) 

t  Time  from  start  of  evaporation/diffusion  (sec) 

x  Space  coordinate  in  slab  geometry  (cm) 

(3  Dimensionless  evaporation/diffusion  coefficient    [Eq.  ( II -8) ] 

|  Integration  variable  (cm) 

T  Dimensionless  time    [Eq.  ( II - 7) ] 
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